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ABSTRACT 

A simple mathematical model was developed to simulate 
the closed portion of the cycle for a quiescent chamber 
compression ignition engine utilizing the assumption of 
perfect gases and the first law of thermodynamics. Various 
input paramentes were used in trend analysis to check the 
model. The output from the computer program was compared to 
test data from a four inch bore, open chamber semi-quiescent 
diesel engine run at the Sloan Test Laboratory. This 
computer model was then modified to simulate the expansion 
stroke of a newly developed, two cylinder cycle 
reciprocating engine. The model was then run to determine 
the optimum point of fuel injection for the new engine. 
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Chapter 1 
INTRODUCTION 



The compression ignition (Cl) engine is thriving in new 
found popularity amongst automobiles, medium-duty and 
heavy-duty freight transport trucks, marine propulsion and 
auxiliary systems, and various other industrial 
applications. The United States Department of Energy (DOE) 
has recently predicted that diesel fuel consumption will 
exceed gasoline consumption in this country by the year 
2000. This is primarily due to the shift to diesel power in 
the automotive and truck freight transport industries to 
take advantage of the high efficiency, high power-to-weight 
ratio engines. The present daily diesel fuel consumption of 

p / 1 \ 

the order of 10 liters is expected to increase by as 
much as 50 percent by the turn of the century. The 
increasing fuel consumption rate, coupled with the constant 
concern of diminishing oil reserves, has prompted renewed 
interest in improving the operating efficiency of the 
conventional compression ignition engine. Any small 
improvement in engine efficiency will obviously result in an 
enormous savings in petroleum. 
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The approaches currently pursued to improve compression 
ignition engine efficiency include increasing the 
compression ratio and the development of the "adiabatic" 
engine. The former involves turbocharging and the 
improvement to piston ring technology. The latter approach 
concentrates the most emphasis on insulating the engine. 
This requires the use of temperature-resistant ceramic 
cylinder liners for combustion cylinders whose gas wall 
temperatures can be of the order of 1200 degrees Kelvin. In 
addition to these two approaches, there are many other 
avenues of research in progress that involve improvements 
that will increase the compression ratio, decrease the heat 
loss from the engine, or increase the combustion efficiency 
through improved combustion chamber design. 

Instead of improving upon the conventional compression 
ignition engine, a new cycle engine design is under 
development. This new design, proposed by Carmichael ^ * , 
consists of a two cylinder cycle which divides the functions 
of a conventional four-stroke diesel cycle into two parts. 
The new engine has one cylinder which compresses the 
incoming air charge and another cylinder which acts as the 
combustion chamber and expansion cylinder. These two 
cylinders are interconnected by a regenerative heat 
exchanger. The regenerator acts as the heart of the new 
design. Through the use of new ceramic materials, the 
regenerator will act as a heat transfer medium by 
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transfering a portion of the heat from the exhaust gases to 
the incoming air charge. The temperature of the incoming 
air will be elevated twice. The first temperature increase 
is due to the compression process in the first cylinder. 
This cylinder, in turn, will transfer its air charge through 
the ceramic matrix of the regenerative heat exchanger, thus 
boosting the temperature for the second time. After passing 
through the regenerator, the incoming air charge will be of 
sufficient temperature to accomodate spontaneous combustion. 
With this high temperature, the high compression ratio of 
the conventional compression ignition engine is not required 
to obtain work from the cycle. Figures 1 and 2 depict the 
pressure vs. volume and temperature vs. volume diagrams for 
the new cycle as compared to a conventional diesel cycle. 
The preliminary design of this new two cylinder cycle engine 
indicates that an improvement to thermal efficiency can be 

achieved over the conventional compression ignition engine. 

An important element of the engine design process is 
the capability to predict, with an acceptable degree of 
accuracy, the energy release during combustion as a function 
of time. This process is extremely complex in that it 
involves the injection and atomization of fuel, the 
evaporation and mixing of the fuel with the air charge, 
followed by the various phases of combustion. The ability 
to accurately predict the heat release rate is vital to the 
engine designer when analyzing a new engine design. 
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This thesis is an attempt to assimilate various diesel 
engine combustion models to produce a simple, yet accurate, 
model to be used in the continuing evaluation of the new two 
cylinder cycle reciprocating engine. The proposed model can 
be utilized on a personal computer system to determine the 
optimum point of fuel injection for the new engine. The 
model has the capability to evaluate two different fuel 
types (i e. iso-octane and propane). 
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Chapter 2 

DESCRIPTION OF THE TWO CYLINDER CYCLE 
RECIPROCATING ENGINE 

The Two Cylinder Cycle Reciprocating Engine consists of 
one compression cylinder and a pair of combustion/expansion 
cylinders, see figures 3 and 4. The compression cylinder 
consists of an intake valve and two exhaust valves, one to 
each expansion cylinder. Each expansion cylinder has its 
own fuel injector. The regenerator cavity contains an 
exhaust valve in addition to the ceramic matrix regenerator. 
The pistons of both the compression cylinder and the 
expansion cylinders are considered to be of simple geometry 
with flat heads. ~ The expansion cylinder piston incorporates 
no unique features to increase turbulence or swirl, thus it 
is similar to a direct injection, quiescent chamber diesel 
engine cylinder. The five valves and three pistons are 
actuated by a camshaft that allows the compression piston to 
operate at twice the speed of an expansion piston. The 
compression cylinder will alternately provide a compressed 
air charge to each expansion cylinder via the regenerator. 
A typical cycle can be illustrated by referring to figures 3 
and 4 . 
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has just closed (it only opened for a very short time just 
before the compression cylinder reached TDC) , valve 3A just 
opened, valves 2A and 3B are closed. (This is the same as 
step 1 except that cylinders A and B are reversed.) 

Step 6 : The compression cylinder is approximately 90° 
ATDC and is in the middle of an air charge induction. 
Cylinder A is still exhausting and cylinder B is in the 
expansion process. Valve 1 is open, valves 2B, 3B, 2A are 
closed and valve 3A is open. (This is the same as step 2 
except that cylinders A and B are reversed.) 

Step 7 ; The compression cylinder has just reached BDC 
and has completed induction of an air charge. Cylinder A is 

exhausting through the regenerator and cylinder B is still 

« 

expanding. Valve 1 just closed, valves 2B, 3B, 2A are 
closed and valve 3A is open. (This is the same as step 3 
except that cylinders A and B are reversed.) 

Step 8 : The compression stroke is approximately 90° 
BTDC and in the middle of compressing an air charge. 
Cylinder A is completing its exhaust stroke and cylinder B 
is ending its expansion stroke. Valve 1 is closed, valve 
2B, 3B and 2A are closed, valve 3A is still open. (This is 
the same as step 4 except that cylinders A and B are 
reversed. ) 
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Step 9 : This is the same as step 1. 



Figures 5 and 6 show the temperature and pressure as a 
function of cylinder volume for a cycle. 

The table below summarizes the sequencing of the valves 
for a complete cycle of an expansion cylinder. 



Tab le -Ij Se quencing Qf valve s 











Valve 
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2A 


1A 


2 B 


2B 


Step 


1: 


X 


X 


X 


X 
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Step 


2: 


0 


X 


X 


X 


0 


Step 


3: 


X 


X 


X 


X 


0 


Step 


4: 


X 


X 


X 


X 


0 


Step 


5: 


X 


X 


0 


X 
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Step 


6: 


0 


X 


0 


X 


X 


Step 


7: 


X 


X 


0 


X 


X 


Step 


8: 


X 


X 


0 


X 


X 


Step 


9: 


X 


X 


X 


X 


0 






where 


X = 


Valve 


closed 








and 


0 = 


Valve 


open 





As can be readily seen, the valve timing sequence is 
rather complex. The timing sequence must be such as to 
allow the impulse transfer of the air charge to occur 
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without possible blow-down to the atmosphere or charging the 
wrong cylinder. A shift of the crank angle must be 
considered to optimize the air charge transfer sequence to 
the on-line expansion cylinder. Thus, the valve timing 
sequence is a critical factor in the correct and efficient 
operation of this new engine design and must be dealt with 
appropriately. 
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Chapter 3 

COMBUSTION AND COMBUSTION MODELING 
(An overview) 

3 . 1 Description of Diesel Engine Combustion 

The diesel engine combustion process is exceedingly 
complex and not very well understood. Combustion in the 
diesel engine is characterized by compression ignition, a 
non-uniform fuel and air distribution in the combustion 
chamber, and a continuous mixing throughout the period in 
which combustion occurs. Due to the initial conditions in 
the chamber when fuel is first injected, the air charge in 
the cylinder is of sufficient temperature and pressure to 
support a chain-reaction. However, combustion in the 
compression ignition engine is governed by the local 
conditions in each part of the charge and not dependent on 
the spread of the flame from one point to another. 
Therefore, the rate of combustion is dependent on the state 
and distribution of the fuel and upon the pressure and 
temperature within the cylinder. 

3.1.1 The Phases of Combustion 
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Ricardo described the diesel engine combustion process 
as taking place in three stages; namely the delay period, a 
period of rapid combustion, followed by burning at a 
controlled rate.*^* Lyn*^* described the burning process in 
three slightly different phases. The first phase is a 
period of rapid combustion which lasts for only three 
degrees crank angle. The second stage is characterized by a 
decreased rate of heat release lasting approximately 40 
degrees crank angle. The third period consists of the fuel 
burning at a very slow rate which may persist through the 
remainder of the expansion stroke. 

A combination of the descriptions of Ricardo and Lyn 
may be more appropriate. The stages of combustion could be 
divided into ignition delay, premixed burning, diffusion 
controlled combustion and the tail of combustion. ^,6 } 
Figure 7 depicts the four stages of combustion in a heat 
release diagram. 

3. 1.1.1 Ignition Delay 

The term ignition delay, or ignition lag, describes the 
time required by the preliminary reactions that occur prior 
to the appearance of flame. The ignition delay is broken 
down into a physical delay and a chemical delay. The 
physical delay period occurs between the beginning of fuel 
injection and the onset of chemical reactions. During this 
period, the fuel is atomized, vaporized, mixed with air and 
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is sometimes 



raised in temperature. This process 

collectively referred to as preparation. The chemical 

delay period immediately follows the physical delay period 
and terminates at inflammation or ignition. This period is 
characterized by chemical reactions starting slowly with 
pre-flame oxidation of the fuel followed by local ignition. 

The ignition delay will vary according to cylinder 
temperature, cylinder pressure, the type of fuel, the 
initial temperature of the fuel, the characteristics of the 
fuel injectors and the turbulence in the cylinder. The 
physical delay is small for light fuels but can become the 
controlling factor for heavy, viscous fuels. The physical 
delay can be significantly reduced by using high injection 
pressures and high turbulence to expedite the breakup of the 
fuel jet. 

Semi-empirical relationships have been developed to 
describe the ignition delay. An estimate for igniton delay 
was developed by Wolfer in 193 8: 

t = 0.44P _1 * 19 exp(4650/T) 

where: t = ignition delay in milliseconds 

P = cylinder pressure in atmospheres 
and T = temperature in degrees K at 
ignition. 
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, ■ 












An estimate by Clarke* in 1970 is quite similar to that by 
Wolf er : 

t = 0-22exp(5500/T)P" 0 * 727 

where: t = ignition delay in seconds 

T * cylinder temperature in degrees K 
and P = cylinder pressure in N/m 2 . 

Still another empirical expression for ignition delay was developed 

{ Ql 

by Spadaccini and Tevelde from experiments for NASA in 1979 with 
diesel fuel in a steady flow facility: 
t = 2.43xl0 _9 P" 2 exp(41560/RT) 

where t = ignition delay in seconds 
P = pressure in atmospheres 
T = mixture temperature in degrees K 
and R = gas constant in atm cm 2 / gmole°K. 

Figure 8 represents the effects of temperature and pressure on ignition 
delay as determined from the estimates by Wolfer. The Spadaccini and 
Tevelde and Clarke relationships yield somewhat similar results. 

When using ignition delay expressions, it must be emphasized that 
differences in engines, fuel properties (especially cetane number) , fuel 
injectors and actual engine temperatures and pressures make the 
calculation rather approximate. These formulas are also very limited by 
their use of bulk temperatures, with no consideration of local 
compositions or temperatures.*^ 
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3. 1.1.2 Premixed Burning 



In the premixed burning stage, flame occurs at one or more 
locations and spreads turbulently. The rate and amount of combustion 
during this stage is directly related to the fuel preparation rate and 
the length of the ignition delay period. Since this stage of combustion 
is one of premixed combustion, little carbon (soot) is produced 
resulting in little radiation heat transfer. However, since the 
combustion rate is so intense, combustion generated noise is controlled 
by this stage of combustion. Figure 9 depicts premixed burning in a 

cylinder. 

3-1. 1.3 Diffusion Controlled Burning 

Once the prepared, or premixed, fuel has burned, the combustion 

process slows down. The combustion rate in this stage will be dominated 

by the rate of local air entrainment. Since the temperature in the 

cylinder is favorable for ignition in this stage, the air/fuel mixing 

process will control the rate of combustion. This preparation of the 

fuel will be governed by the turbulence and swirl in the cylinder. 
{4} 

Lyn estimated that approximately 40 percent of the heat release frcm 
the combustion of fuel comes frcm this stage. Figure 10 shows the 
diffusion burning process in a cylinder. 

3. 1.1.4 Combustion Tail Stage 

This last stage of combustion is characterized by the cylinder 
pressure and temperature falling as the expansion process continues. 
The rate of combustion tails off due to the chemical kinetic effects as 
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the chemical reaction rate slows. In this stage, the reaction rate will 
become the controlling factor instead of the air/ fuel mixing process. 
This stage is also characterized by diffusion combustion with a high 
production and combustion of soot particles with a resultant high rate 
of radiation heat transfer. This last stage of combustion can proceed 
through the completion of the expansion stroke and can contribute upto 
20 percent of the total heat release. ^ Figure 11 represents a typical 
heat release rate diagram showing the four stages of combustion. 

3.2 Combustion Modeling 

The combustion process is often considered the most important 

aspect of an internal combustion engine, but, at the same time, the 

least understood and most complex. A mathematical model depicting 

combustion would require good models of the fuel system to include the 

injection/fuel pump, the injector nozzles, and fuel lines. 

Additionally, models of fuel atomization, vaporization, fuel/air mixing, 

cylinder air motion, chemical kinetics and pre-mixed and diffusion 

mixing would be required. A model as comprehensive as this has yet to 

{ 12 } 

be deveoped. Spaulding states that this type of "combustion 

modeling is impossible." He justifies this by pointing out that the 
number of governing restraints and rules outnumber the degrees of 
freedom and, in addition, the requirements of low cost, speed and 
accuracy must also be met. Since the complexity of the real combustion 
process is so overwhelming, substantial simplifying assumptions must be 
made to obtain solutions. 
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3.2.1 Types of Models and Uses 



113 } 

Bracco categorized combustion models into three categories 
based on their uses in examining different engine problems. The 
categories are the zero-dimensional (or thermodynamic) model, the 
quasi-dimensional (or entrainment) model, and the multi-dimensional (or 
detailed) model - 

3. 2. 1.1 Zero-dimensional Model * 11 ^ 

The zero-dimensional model is structured around a thermodynamic 
analysis of the engine cylinder contents during the cycle. The 
assumptions include one-dimensional flow, isentropic adiabatic flow 
through nozzles simulating flow past valves, and unburned mixtures as 
mixtures of air, fuel vapor and residual gases. Specific heats of the 
gas mixture are modeled using polynomial functions of temperature. 
Compression is assumed to be adiabatic. Combustion assumes 
thermochemical equilibrium and progressive burning via mass elanents. 
The expansion process assumes thermochemical equilibrium. 

Heat transfer is modeled using correlations between the Nusselt, 
Prandtl- and Reynolds numbers from heat transfer in steady turbulent 
flow over flat plates and pipes. These relationships are in the form 
of: 

Nu = aRe^Pr c 

where a, b, and c are obtained from experimental data for a specific 
engine. 
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The combustion process is generally modeled from an apparent heat 
release or an experimentally obtained fuel burning rate. One of the 
most widely used correlations is based on the Wiebe Function. In this 
function, the fuel burned is expressed as a fraction of the total fuel 
injected. ^ 

FB = 1 - exp[-K 2 (t) (K 1 +1) ] 

where FB = fraction of fuel burned/total 
injected 

t = time from ignition 

= shape factor for combustion curve 
~ combustion efficiency coefficient. 

Another typical function form is the cosine function: 

X(6) = (1/2) (1 - cos Tit (0 - e^/A©^} 

where X(0) = mass fraction burned at crank 
angle 0 

0 Q = crank angle at the start of 
combustion 

and A©^ = burn duration. 

There are numerous other combustion models that utilize various 
heat release patterns. Sane replace the heat release curve with two 
straight lines. In this type of combustion model, one line simulates 
the rapid combustion of the bulk of the injected fuel and the other line 
represents the slower combustion phase further down the expansion 
stroke. 
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An empirical model developed by Whitehouse and Way*^ is based on 
elementary combustion principles. Fuel is assumed to be prepared for 
combustion as a result of fuel-air mixing. The reaction rate calculates 
the burn rate in the premixed stage of combustion. The preparation rate 
becomes governing during the diffusion burning phase as the fuel is 
assumed to burn as rapidly as it is prepared. (The Whitehouse and Way 
model will be dealt with in detail in a later chapter.) 

In general, thermodynamic combustion models are useful when 

performing a design trade off or comparison analysis to evaluate the 

effects of change in engine design and operation. Since, however, the 

details of the combustion process are an input to the model, the results 

can only indicate what will transpire if the engine burns in the 

specified manner. These models cannot address the feasibility of the 

engine operating in the prescribed manner because the details of the 

(15*} 

burning process are not linked to the engine design and operation. 

3.2. 1.2 Ouasi-dimensional Model *-*--*-* 

Quasi-dimensional models are also structured around a thermodynamic 
analysis of. the engine cylinder during the cycle. Many of the same 
assumptions are utilized to describe the various portions of the process 
as are used in the thermodynamic model. The combustion process, on the 
other hand, is based on more fundamental physical quantities such as 
turbulent intensity, turbulent mixing, jet characteristics in jet mixing 
and the kinetics of the fuel-oxidation process. 
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The quasi-dimensional models can be utilized for the same purposes 
as the zero-dimensional models except that they can new be used where 
changes in the combustion process can be a dominant factor. The major 
drawback of the quasi-dimensional model is its inability to examine, in 
detail, the interaction between fluid flow and engine geometry. 

3.2. 1.3 Multi-dimensional Model 

In a multi-dimensional model, the governing partial-differential 
equations describing conservation of mass, momentum, energy and species, 
and the sub-models describing turbulence, chemical kinetics, and etc. 
are numerically solved subject to boundary conditions and other 
restraints. These models have the potential for examining the 
interaction between fluid flow and engine geometry that is lacking in 
the quasi-dimensional model- The detailed model will predict engine 
performance and emission characteristics from the first principles with 
virtually no empirical relationships. Unfortunately, solving the 
relevant conservation equations in three-dimensional, time dependent 
formulation, coupled with the state equations and sub-models leads to a 
computer program that will tax even the most capable computer system. 
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Chapter 4 

HEE TWO CYLINDER CYCLE C08COSTION MCDEL 



Since the two cylinder cycle reciprocating engine is a totally new 
concept, combustion modeling can be even more difficult than for a 
compression ignition engine. However , the approach taken models the 
expansion cylinder of the new cycle after a diesel engine cylinder. The 
beginning of the expansion stroke will simulate a diesel engine with its 
piston at TDC with a charge of air. For this initial combustion model, 
the air will be assumed to be contained within the cylinder, at 
pressure, with no additional air added after expansion, as in the actual 
new engine cycle. 

4.1 Assumptions 

The assumptions for this single zone combustion model are 
essentially those previously mentioned for the thermodynamic type of 
models. 

a. The First Law of Thermodynamics is used to establish an energy 
balance to determine the temperature at the end of each step. 

b. The working fluid is treated as an ideal gas. 

c. The system contents are homogeneous and of uniform temperature 
and pressure. 
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d. The changes in gas properties due to the rate of change of the 
gas composition are considered to be negligible. 

e. Combustion is treated as a reversible heat release process. 

f. Combustion products are formed in the proportions according to 
the law of perfect combustion. 

g. No dissociation of the products of combustion occurs. 

h. CWLy four gases are considered to be present and are varied as 
required for perfect combustion. 

i. The incoming air charge is assumed to be pure air plus a 
fraction of the residual gases remaining in the cylinder. 

4.2 Thermodynamics of Internal Combustion Engines 

4.2.1 Idea l , £as {16} 

The assumed thermally ideal gas obeys the state equation 
pV = MRT 

where p = pressure 
V = volume 
M = number of moles 
R = universal gas constant 
and T = temperature. 

The specific gas constant, R, can be written in terms of R and m , 

w 

the molecular weight of the gas. 

R= 
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If the mass of the gas, rrFttr^, then the state equation can be written 
as: 

pV = mKT. 

The specific internal energy for an ideal gas can be represented as 
a function of temperature: 
u = f (T) 

where u = specific internal energy 
and f (T) = function of temperature dependent on 
the gas. 

If the function f(T) is expressed in the form of a limited power series, 
then' 171 



Thus, following the same procedures as for the internal energy, 




where a^ to a^ are constants which vary 
depending on the gas 

and u Q •= internal energy at absolute zero. 



The specific heat at constant volume can be defined as: 



C v = (dq/dT) v =(du/dT) v . 




rvi| 

The specific enthalpy, h, for an ideal gas is given by: 
h = u + RT. 

It follows that: {17 ' 18} 



n=£ 



h = h(T) = u Q + / + RT 
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At absolute zero f T=0: 
h=h o =u o* 

Therefore, for a perfect gas, the internal energy varies linearly with 
temperature as: 

h = h Q +CyT + RT. 

The specific heat at constant pressure, C pf is defined by: 

C p = (dq/dT) p = (dh/dT) p . 

For a perfect gas: 

C p - Cv + 

New, enthalpy can be expressed by: 
h - h o + c p T - 

For thermodynamic processes with gases of constant composition and 
specific heats undergoing state changes; 

h o ■ u o ■ °- 

Then, 

u - C V T; 

h - CpT; 

h - u - (Cp - CyJT = RT; 
and C p - = R. 

Gas data are often given in terms of enthalpy vice internal energy. 
The conventional form is: 

h(T)/RT = (h - h )/RT 

= a 1 + ^T + a 3 T 2 + a 4 T 3 + a 5 T 4 . 
and the internal energy is expressed as: 

u(T)/RT = (a.j-1) + c^T + a 3 T 2 + a 4 T 3 + a $ T 4 . 
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The values for the polynomial coefficients, a^ to a^ are provided in 

Table 2. Other formulations for the calculation of enthalpy and 

{ S 26 27 } 

specific heat are available in the literature. ' ' 



Table 2: Polynomi al Coefficients 

Range: 500 - 3000 Degrees Kelvin 







^2 


% 


*4 


% 


C0 2 


3.0959 


2.73114E-03 


-7 . 88542E-07 


8.66002E-11 


0.0 


* 2 ° 


3.74292 


5.65590E-04 


4.95240E-0 8 


-1.81802E-11 


0.0 


°2 


3.25304 


6.52350E-04 


-1.49524E-07 


1.53 897 E-ll 


0.0 


*2 


3.34435 


2.94260E-04 


1.95300E-09 


-6.57470E-12 


0.0 


C 8 H 18 


-0.71993 


4.6426E-02 


-1.68385E-05 


-2.67009E-09 


0.0 


C-jHp 


1.13711 


1.45532E-02 


-2.95876E-06 


0.0 


0.0 



4.2.2 Pro perties of Gas Mixtures 1X0 J 

Mixtures of gases obey the following. 

a. The gas mixture as a whole obeys the equation of state, 
pV = MRT, where M is the total number of moles of all species. 

b. The total pressure of the mixture is equal to the sum of the 
pressures which the individual ccmponents/species exert. 

c. The internal energy, enthalpy and entropy of the mixture equals 
the sum of the internal energies, enthalpies and entropies which each 
individual component/species would have if it separately occupied the 
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entire volume of the mixture at the same temperature 



Thus, for mixtures of ideal gases the mole fraction is given by: 



x i = Mj/M 



where = moles of a specie 
and M = total number of moles. 



Then, 



Xx i = l.o. 



Enthalpy is given by: 

H = XM i h i = M Z x i h i . 
Internal energy is given by: 

U = Xm.^ = M £ x^u i . 
Specific Heats are given by: 

c p = S x i c pi 
c v " = x i c vi- 

4.2.3 The First Law of Thermodynamics 



{17} 



The emphasis of this model is the closed portion of the cycle. 
Therefore, the First Law of Thermodynamic for closed systems is simply: 
dQ - dW = dU' 

where dQ = heat energy transfer 
dW = work energy transfer 
dU' = change in internal energy. 

The internal energy is defined by: 

U' = U + KE + PE 

where U = the intrinsic internal energy 
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KE = kinetic energy 
PE = potential energy. 



For a closed system, we can assume that PE = KE = 0. 
dQ - dW =dU 

where U = M I x^u^ 

M = total number of moles 

x^ = mole fraction of gas i 

u^ = specific internal energy of gas 

For non- reacting closed systems, we can write: 
dQ - dW = dU 

where dW = pdV = ( I x^p) dV 
and dU = Md( I x^u^) . 

For a reacting closed system, we can expand this to: 

dQ - pdV = (U - U or ) + U (T) - O r (T) 

where (U - n ) = ^ u 

op or o 

1 U Q = heat of reaction 
U (T) = energy of products as a 
function of time 

U r (T) = energy of reactants as a 
function of temperature. 

4.3 Heat Transfer from the gas to the Cylinder 



Therefore, 



To be able to balance the energy in a real system, the heat 
transfer f ran the combustion gas to the walls of the cylinder must be 
considered. TWo basic equations are generally accepted for use in cycle 



calculations. These are the correlations developed by Annand and 

Woschni. The relationship by Woschni {19} is based upon a forced 
convection model. 

q/A=C 3 d~° * 2 p° * * 053 (C 1 V p + (C 2 (p-p o ) VT ' /p ' V ' ) 0 * 8 (T g -T w ) 

where C^, C,, and C 2 = constants 
A = area 

D = cylinder bore 
p = pressure 

Tg = mean gas temperature 
T w = wall temperature 
V p = piston velocity 
p Q = motoring pressure 
p' = trapped pressure 
V* = trapped volume 
T* = trapped temperature. 

Although Woschni* s expression is readily accepted, it does not 
separately distinguish between convection and radiation. 

The Annand equation is also largely based on turbulent convection. 
Unlike the Woschni correlation, Annand claims that the Reynolds number 
is the major parameter affecting convection. Convection is the first 
term in his equation. The second term in Annand *s equation is a 
radiation term assuming grey body radiation. Thus:* 29 * 
q/A = a(k/D) (Re) b (T g -^) + c(Tg 4 —, IW 4 ) 
where q = heat transfer rate 
A = area 
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a,b,c = constants 
k = thermal conductivity 
D = bore 



Re = Reynolds Number = pVjrD/u 
p = density 
Vp = piston velocity 
P = viscosity 

Tg = temperature of gas (mean) 

Tw = temperature of wall. 

The range of values for Annand's constants are: 
for a four stroke engine: 
a = 0.26 
b = 0.75 ± 0.15 

c = 3.88 ± 1.39 x 10" 8 J/sm 2 K 4 
for a two stroke engine: 
a = 0.26 
b = 0.64 ± 0.10 
c = 3.03 ± 1.06 x 10" 8 

Since Annand's equation separates the convective term 
radiation term, it is believed that the Annand correlation 
suited to the new cycle calculations. 

4.4 The Combustion Model 



f ran the 
is better 
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In the process or heat release from combustion, both physical and 
chemical effects are involved. Liquid fuel injected into an engine 
must be heated, vaporized, and mixed with oxygen in the preparation 
process prior to combustion. Once the fuel is prepared, it may then 
burn at a rate controlled by chemical kinetics. It has been 
demonstrated that the time required for combustion of the prepared fuel 
is negligible as compared to the preparation time. 

At the beginning of the burning period, chemical kinetics are 
important due to the low temperatures. When fuel is first injected into 
a cylinder of a diesel engine, the temperature is generally such that 
rapid burning will not occur. Additionally, the heat transferred to the 
incoming fuel causes the temperature to drop in the cylinder. As the 
temperature rises in the cylinder, the combustion rate rises, thus 
increasing the temperature. The heat release rate continues to rise 
until the lack of prepared fuel becomes the controlling factor. When 
the excess prepared fuel is depleted, combustion will proceed at the 
rate of fuel preparation. Figure 12 represents the effects of 
preparation rate and reaction rate in premixed burning as a function of 
crank angle. 

4.4.1 Preparation of Fuel 

After injection, the fuel is physically prepared for combustion. 
As mentioned before, this process involves the atomization, vaporization 
and mixing of the fuel with air. The rate of preparation can be assumed 
to be proportional to the total surface area of the fuel spray droplets. 
If all the droplets are assumed to be of identical size, then it 
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follows: < 7 ' 14 ' 21 > 

M^ = np /r D 0 2 /6 
^ = np n D 2 

where M^ = Mass of fuel injected 
= Mass of fuel unburned 
n = number of fuel droplets 
p = fuel droplet density 
D 0 = Initial droplet diameter 
D = Droplet diameter. 

The total area 

Area = n^D 2 = n n (6M u /np:r J 273 
Area = ^/apD^) 173 ^/p*) 273 
Area = 6 M^ 7 3 M^ 27 3 / pD Q . 

Assuming that the density, p, and initial diameter, D Q , are constant, 
then the 

4 

Area a M.^ 73 M 273 . 

" l u 

Allowing for the effect of oxygen availability on the mixing of the 
fuel, the preparation rate, PR, can be written as: 

PR = KM i 1_x M u x P 02 m 

where x = empirical constant 
m = empirical constant 
P 02 = partial pressure of oxygen 
K = constant. 

The constant K is a function of the characteristics of fuel injection, 
air movement and combustion chamber shape. Typical values for four 
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. ( 14 ) 

stroke engine are: 

K = 0.008 - 0.020 

x = 2/3 

m = 0.4. 

4.4.2 Reaction of Fuel 

Since diesel fuel is not a pure substance, it is impossible to 
ascertain the exact chemical equations involved since the actual 
compounds in the fuel are unknown. The temperatures that are available 
f ran experiments are only average cylinder temperatures. With these 
approximations/estimations, the equations for reaction rate are highly 
empirical. The degree of approximation involved may be justified due to 
the short time period during which chemical kinetics is of importance. 
Also, the total fuel that is burned is equal to the amount of fuel that 
is prepared. The reaction rate equation that was proposed by Whitehouse 
and Way * 7 ' 14 ' 21 * is based on the Arrhenius equation. 

R = (K'P^j/ftvr” ) J( PR-R)dx exp(-act/T) 

where R = reaction rate per degree crank angle 
K' = empirical constant 
act = empirical constant 
Pq 2 = partial pressure of oxygen 
PR = preparation rate 
N = engine speed in rpm 
T = cylinder temperature. 

The effect of the ignition delay period is incorporated in the Arrhenius 
type expression exp(-act/T). Typical values of K' and act are : 
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act = 1.4 x 10 4 

K' = 1.2 x 10^® for two stroke engines 
K' = 65 x 10^® for four stroke engines 
4.5 Verification of the Model 

The model was converted top computer code using TRS-80 Model III 
Disk Basic. The program listing is presented in Appendix B. 

In an effort to set the empirical coefficients, the average value 
was used for all coefficients that had a range of values for four stroke 
engines. The program was run and compared to the data obtained by 

Remley in actual engine testing in the Sloan Automotive Laboratory. 

Figure 13 represents the pressure versus volume curve for the model and 
for the engine run by Remley. Appendix A provides specifications of the 
test engine. 
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Chapter 5 

SELECTION CF FUEL INJECTION POINT 



In order to obtain the maximum work and highest efficiency frcm the 
new two cylinder cycle, the time of fuel injection should be optimized. 
To obtain this optimum, a number of cycles were run on the computer. 

5.1 Selection of the Model Coefficients 

The model was run assuming the expansion cylinder at ICC with an 
air charge at a temperature and pressure of 1090°K and 10 atmospheres 
while the engine speed of 850 rpm and air/fuel ratio were held constant. 
The selected fuel was CgH^g (iso-octane) with a lower heating value of 

4.2 x 10 7 joules/kilogram and a residual air fraction of 0.05. 



The model was run several times to obtain a value of K in the 
equation: 

PR = KM i (1-x) M u x P 02 . 

The values of x and m were held constant at 2/3 and 0.04, respectively, 
as the values used for four-stroke diesel engines. When searching for a 
value of K, a diffusion combustion period of 70 - 120 degrees of crank 
angle was sought. This was found through several iterations to occur at 
a value of K = 0.012. 
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The values of constants for the reaction rate equation: 

R = [K'P 02 /NT°* 5 ]exp(-act/T) J(PR-R)dx, 
were selected as the values for four-stroke diesel engines. 

With this input data and selection of constants, the model yields a 
heat release rate curve which closely resembles that described by 
Ricardo, Lyn and Whitehouse et al, see figure 14. The premixed burning 
phase yields approximately 45 percent of the heat release, the diffusion 
controlled burning phase yields approximately 45 percent of the heat 
release with the tail of combustion providing the remaining 10 percent. 

5.2 Optimizing Fuel Injection 

Intuitively, the maximum work and highest efficiency would be 
expected with fuel injection and combustion occur ing at TDC, or 
immediately thereafter. This, however, does not appear to be the case 
when the data is evaluated. See figures 15 through 21. While the fuel 
injection is varied from 180 (TDC) to 205 degrees crank angle, with an 
injection period of 20 degrees, the thermal efficiency rises. For fuel 
injection occurring from 180 to 195 degrees, the temperature at 360 
degrees (BDC) is not sufficient to heat the regenerator matrix to a 
temperature which will pre-heat the incoming air charge to 1090 degrees 
Kelvin as specified by the input data. For fuel injection occuring at 
205 degrees, and later, incomplete combustion will result. 

From this approach, the optimum point of fuel injection occurs at 
200 degrees crank angle for a fuel injection period of 20 degrees. 
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When the fuel injection period is reduced to 10 degrees, a similar 
pattern is observed. A fuel injection point on, or before, 200 degrees 
results in the cylinder gas temperature dropping too low to support 
sufficient air charge pre-heat. Fuel injection on, or after, 210 
degrees results in incomplete combustion. See figures 22 through 24. 
In this case, the optimum point of fuel injection occurs at 205 degrees. 
The thermal efficiency for this case is higher that the case of a 20 
degree injection period. Also, the specific fuel consumption is lower 
in the case of 10 degree injection as compared to 20 degree injection. 

Through a similar analysis, the case of an air/ fuel ratio of 25 
yields an optimum fuel injection point of 195 degrees crank angle for a 
period of 20 degrees. For this air/fuel ratio, the value of K in the 
preparation rate equation was selected as 0.018 to achieve a similar 
heat release rate pattern. 
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Chapter 6 

COMMENTS AND RECOMMENDATIONS 



As can be readily seen, the output from this type of thermodynamic 
model is dependent on the value of the empirical coefficients. The 
characteristics of the heat release rate curve will shift as a function 
of air/fuel ratio, temperature, pressure and engine speed. Therefore, 
only a comparison of results from a defined heat release should be used 
for qualitative comparison analysis. 

Since this computer model was written for a personal computer, the 
time required for one run is excessively long for a detailed comparison 
analysis. The run time for one run with five degree increments is 
approximately 1 hour and 20 minutes. A motoring analysis (no 
combustion) requires approximately 15 minutes. The amount of time 
required in the combustion iteration process is the difference between 
the two. These times were obtained when running the program with no 
remark statements and elimination of all unnecessary spaces in the 
program. Undoubtedly, the efficiency of the program can be somewhat 
increased by utilizing some clever programming techniques. However, the 
use of a small computer strictly dedicated to a comparison analysis with 
crank angle increments of one or two degrees can occupy the machine for 
an inordinate period of time. 



- 43 - 



6.1 Recommendations 



This single zone model allows for cycle studies. However, a 
problem that must be explored is the formation of soot and gaseous 
pollutants. This can be accomplished by expanding the model to a two or 
four zone model. *23 '24} D Ur j. n g this model expansion, the effect of 
chemical kinetics should be further examined to display a more realistic 
combustion process. The values of the coefficients for the polynomial 
expression of enthalpy for the other products of combustion are readily 
available . * 17 

The effects of heat transfer from the system may be more 
appropriately modeled by the use of the widely accepted Woschni 
correlations. The use of Annand's correlation, however, does allow 
for the separation of convection and radiation. 

The effects of mixing of the air charge with the fuel must be 
further explored to determine the effects on combustion intensity and 
efficiency. 

The use of a larger computer system would be most beneficial in a 
comparison analysis. Single runs can be easily done on a personal 
computer system, however, many runs using small crank angle increments 
are best, although more costly, performed on a main frame system capable 
of performing numerous simultaneous calculations. 

Lastly, to obtain realistic coefficients for the empirical 
constants in the preparation and reaction rate equations, experiments 
using a rapid compression machine are considered appropriate. This 
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would provide for realistic data with minimum cost. 
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Figure 1: Pressure vs. Volume for New Cycle vs. Oiesel Engine 
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Figure 2: Temperature vs. Volume for New Cvcle vs. Diesel 
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Figure 5: Temperature vs. Volume for Mew Cvcle 
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Figure 6: Pressure vs. Volume for New Cvcle 
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Figure 7: The Four Phase? of Combustion 
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Figure 8: Effects of Temperature and Pressure on Ignition Delay 
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Figure 10: Direct-Ini ection Compression Ignition Engine 

Combustion during the MIXING CONTROLLED Phase 
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Figure 14: Heat Release Rate Obtained from Computer 
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DIESEL ENGINE COMBUSTION CYCLE 



RUN: CARMICHAEL ENGINE 
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DIESEL ENGINE COMBUSTION CYCLE 



RUN: CARMICHAEL ENGINE 

INPUT data: 

CYLINDER BORE = .3725 METERS 

STROKE = .3725 METERS 

CONNECTING ROD LENGTH = 

ENGINE SPEED = 850 RPM 

ENGINE COMPRESSION RATIO = 5 

AIR / FUEL RATIO = 30 

TRAPPED PRESSURE = 1.01325E-06 

TRAPPED TEMPERATURE = 1090 DEG 

RESIDUAL AIR FRACTION = .05 

FUE1I SELECTED FOR THIS ANALYSIS 
WITH A LOWER HEATING VALUE 
STOICHIOMETRIC AIR / FUEL RATIO = 15 



,745 METERS 



N/MT2 

KELVIN 



= C8H1S (ISO OCTANE) 
= -4.2E+07 JOULES/KG 
.1151 



FUEL / AIR EQUIVALENCE (PHI) = 



,503838 
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DIESEL ENGINE COMBUSTION CYCLE 



RUN: CARMICHAEL ENGINE 

INPUT DATA: 

CYLINDER BORE = .3725 METERS 

STROKE = .3725 METERS 

CONNECTING ROD LENGTH = .745 METERS 

ENGINE SPEED = 850 RPM 

ENGINE COMPRESSION RATIO = 5 

AIR / FUEL RATIO = 30 

TRAPPED PRESSURE = 1 . 01325E+06 N / M 1 2 

TRAPPED TEMPERATURE = 1090 DEG KELVIN 

RESIDUAL AIR FRACTION = .05 

FUEL SELECTED FOR THIS ANALYSIS = C3H18 
WITH A LOWER HEATING VALUE = -4.2E 
STOICHIOMETRIC AIR / FUEL RATIO = 15.1151 

FUEL / AIR EQUIVALENCE (PHI) = .503836 



(ISO OCTANE) 
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DIESEL ENGINE COMBUSTION CYCLE 





RUN: CARMICHAEL ENGINE 
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DIESEL ENGINE COMBUSTION CYCLE 



RUN: CARMICHAEL ENGINE 



,745 METERS 



INPUT data: 

CYLINDER BORE = .3725 METERS 

STROKE = .3725 METERS 

CONNECTING ROD LENGTH = 

ENGINE SPEED = 850 RPM 

ENGINE COMPRESSION RATIO = 5 

AIR / FUEL RATIO = 30 

TRAPPED PRESSURE = 1.01325E+06 

TRAPPED TEMPERATURE = 1090 DEG 

RESIDUAL AIR FRACTION = .05 

FUEL SELECTED FOR THIS ANALYSIS 
WITH A LOWER HEATING VALUE 
STOICHIOMETRIC AIR / FUEL RATIO = 15 



N/MT2 

KELVIN 

= C8H18 
= -4.2E 
.1151 



'UEL / AIR EQUIVALENCE (PHI> = .503836 



(ISO OCTANE) 
07 JOULES/KG 
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DIESEL ENGINE C0M3US7I0N CYCLE 



RUN: CARMICHAEL ENGINE 

INPUT DATA*. 

CYLINDER EGRE = .3725 METERS 

STROKE = . 3725 METERS 



CONNECTING ROD LENGTH = .745 METERS 

ENGINE SPEED = 350 RPM 

ENGINE COMPRESSION RATIG = 5 

AIR / FUEL RATIO = 30 

TRAPPED PRESSURE = 1.01325E-06 N/Mt2 

TRAPPED TEMPERATURE = 1090 DEG KELVIN 

RESIDUAL AIR FRACTION = .05 

FUEL SELECTED FOR THIS ANALYSIS = C8H18 
WITH A LOWER HEATING VALUE = -4.2E 
STOICHIOMETRIC AIR / FUEL RATIO = 15.1151 

FUEL / AIR EQUIVALENCE (PHI) = .503836 



(ISO OCTANE) 
07 JOULES/KG 




DIESEL ENGINE COMBUSTION CYCLE 



INPUT DATA: 



RUN: CARMICHAEL ENGINE 

CYLINDER BORE = .3725 METERS 

STROKE = .3725 METERS 

CONNECTING ROD LENGTH = .7^5 

ENGINE SPEED = S50 RPM 

ENGINE COMPRESSION RATIO = 5 

AIR / FUEL RATIO = 30 

TRAPPED PRESSURE = 1.01325E 



METERS 



'RAPPED 



iMPERATURE 



06 

1090 DEG 



N/MT2 

KELVIN 



RESIDUAL AIR FRACTION = .05 

FUEL SELECTED FOR THIS ANALYSIS 
WITH A LOWER HEATING VALUE 
STOICHIOMETRIC AIR / FUEL RATIO 
FUEL / AIR EQUIVALENCE (PHI) = 



C3H18 
-4.2E 
15.1151 
503836 



(ISO OCTANE) 
07 JOULES/KG 
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Figure 22 
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ENGINE COMBUSTION CYCLE 



DIESEL 



RUN 

INPUT data: 

CYLINDER 
STROKE = 



CARMICHAEL ENG INI 



BORE = .3725 

.3725 METERS 
CONNECTING ROD LENGTH = 
ENGINE SPEED = 850 RPM 



METERS 



.745 METERS 



RATIO = ‘ 

30 

1 . 0 1 325E J 



ENGINE COMPRESSION 
AIR / FUEL RATIO = 

TRAPPED PRESSURE = 1.01325E-06 

TRAPPED TEMPERATURE = 1090 DEG 

RESIDUAL AIR FRACTION = .05 

FUEL 



SELECTED FOR THIS ANALYSIS = 
WITH A LOWER HEATING VALUE = 
STOICHIOMETRIC AIR / FUEL RATIO = 15. 

FUEL / AIR EQUIVALENCE (PHI) = .50383 



/ M f 2 

ELVI N 

CSH1S 
-4 . 2E 
1151 
6 



(ISO OCTANE) 
07 JOULES/KG 
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DIESEL ENGINE COMBUSTION CYCLE 



RUN: CARMICHAEL ENGINE 

INPUT DATA.’ 

CYLINDER BORE = .3725 METERS 

STROKE = .3725 METERS 

CONNECTING ROD LENGTH = .745 METERS 

ENGINE SPEED = 850 RPM 

ENGINE COMPRESSION RATIO = 5 

AIR / FUEL RATIO = 30 

TRAPPED PRESSURE = 1.01325E + 06 N / M f 2 

TRAPPED TEMPERATURE = 1090 DEG KELVIN 

RESIDUAL AIR FRACTION = .05 
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Appe ndix h 



Specifications of Test Remley Engine 



{ 22 } 



Type of Engine 

Bore 

Stroke 

Cylinder Displacement 
Connecting Rod Length 
Compression Ratio 
Number of Compression Rings 
Number of Oil Rings 
Number of Inlet Valves 
Number of Exhaust Valves 
Valve Diameter 
Valve Lift 

Inlet Valve Timing open/close 
Exhaust Valve Timing open/close 
Diameter of Intake Manifold Pipe 
Diameter of Exhaust Manifold Pipe 



Four Stroke 

4.0 inches 
2.5 inches 

31.41 cubic inches 
6.25 inches 
14.3 : 1 
2 
1 
2 
2 

1.286 inches 
0.280 inches 
1 5 °BTDC/ 50 °ABDC 
50°BBDC/15°ATCC 

2.00 inches 
1.60 inches 
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Data Collected from Test Run 



Engine Speed in RPM 
Inlet Pressure 
Exhaust Pressure 
Inlet Air Temperature 
Air to Fuel Ratio 
IMEP 

Start of Injection 

Ignition Delay 

Period of Fuel Injection 



{22J 

-1450 

13.5 inches Hg gage 

13.4 inches Hg gage 

186° F 

25.38 

88.1 psi 

12.5° BTCC 

5.5° 

17.5° 
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Appendix B 

Computer Model 



The computer program is written in TRS-80 Model III Disk Basic and 
consists of a main program and nine subroutines. The program listing 
has numerous remarks statements inserted to make the algorithm and 
computer code easier to understand. Since the program takes a 
considerable length of time to run, it is recommended that the remark 
statements be deleted before running. Samples of output are presented 
in figures 15 to 24. 
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This program is written in TRS-80 Model III Disk Basic. 

Remove all remarks spaces before running to speed up run time 
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160 INPUT“ENTER CRANK ANGLE FOR EXHAUST VALVE 0PEN 4, ;AEV0 
170 INPUT-ENTER CRANK ANGLE FOR FUEL INJECTION -;AIJECT 



ISO INPUT-ENTER PERIOD OF FUEL INJECTION ( DEGREES) I JECT 
185 INPUT-ENTER CRANK ANGLE INCREMENTS FOR THIS RUN" ; ADELT 

190 INPUT-SELECT FUEL: <1> FOR C8H18 (OCTANE) (2) FOR C3H8 (PROPANE) 

195 'Load in fuel data for selected fuel. 

200 IF NN= 1 THEN GOSUB 4000 ELSE GOSUB 4100 



z 

z 

• ft 























Q 






















LQ 






















O 






















Z 






















Ui 






















Z 






















z 






















o 






















u 






















z 






















o 






















M 






















H 












c 










cn 












0 










p> 






















CQ 












z 










Z 












0 










O 












p 










O 












z 






















£ 






















0 
















• 






a 










O 






0 
















ui 






0 






z 










*-> 






P 






u 










M 






(0 






0 










< 












V 










• ft 






0 






u 










z 






a 






0 
















0 






a 










H 






vi 
















< 






p 






c 
















p 






0 










H 






CO 
















a 












V 










< 






V 






0 










H 






c 






0 










cn 






0 






0 






















z 










z 






0 
















o 






£ 
















H 






P 






0 










H 






•H 






P 










U 


o 


« 


0 






Vl 










Ui 


o 




> 


















(N 










Vi 










z 


Tf 




p 






0 










M 






0 


















0 




*0 






0 










Z 


C 




c 






0 










LU 






•H 






rH 










D 


-P 




•H 






0 










Ui 


3 










£ 










■ 


0 




a 
















ft 


P 










Vi 










% 


A 




0 






0 










H 


P 




0 
















Z 


0 




Z 






P 










w 






0 






0 










OC 


Z 




•—4 






Z 










CU 


-p 




P 






£ 




■ 






Z 


•H 




u 






P 




0 








3 










C 




Z 






z 






0 










0 






Ui 


10 




0 






•o 




-o 






X 


Z 










c 










H 


C 




o 






0 




z 








0 




o 










p 






< 


-P 




o 






*H 




a 






X 


10 




n 






z 




c 






Cu 


c 










a 




•H 




cc 


Z 


0 




0 














< 


< 


0 


o 


C 


O 




0 


o 


Z 


o 


CQ 


it 




o 


•h O 




z 


o 


P 


o 


CU 


H 


c 


CM 


Z O 




0 


rH 


0 


CM 


* 


O 


*H 


TP 


P 


in 




•H 


m 




in 


rH 


Ui 






0 






P 




z 




OU 


*-> 


V 


CO 


u 


OQ 




u 


CQ 


c 


OQ 


II 


W 


0 


ZD 


Z 


D 


> 


rH 


p> 


•H 


ZD 


cn 


< 


0 


cn 


p cn 


II 


0 


cn 


p 


cn 


ID 




J 


o cn o 


r-l O 


o 


a 


o 


CC 


Ui 


X 


o 


\ 


a 


> 


\ 


o 


ft 


o 


CU 


W 


in 


o 


in 


o 


in 




o 


in 


o 


o in 


o 


rH 


rH 


CM 


CM 


CM 


CO 


CO 




in 


in 


(N 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 



• 0 

C *0 
p -H 









P 


X 


U 
















• 






0 


0 
















0 




V 


♦ri 


a 
















rH 




0 




0 




c 












0 








> 




0 












£ 




cn 


C 




C 


cn 












CO 




c 


0 


p 


0 


0 












X 




-rH 


XJ 


0 


CT> 


p 


rH 










\ 




c 


u 


-p 


>S -p 


0 










■-> 




c 


0 


0 


X 


•H 


P 














•H 


O 3 


o 


z 


u» 










0 
























z 




0 


V 


V 


V 


V 


V 














z 


0 


0 


0 


0 


0 








o 


cn 






















x 


x 




4j 


0 


0 


0 


0 


0 








ft 


\ 




0 


0 


0 


0 


0 


0 








u« 


*-> 






rH 


rH 


rH 


rH 










z 






0 


0 


0 


0 


0 


0 








a 


£ 




0 


z 


z 


z 


z 


z 










0 




•H 


N 


s 


ft 


ft 


ft 








H 


p 




CJ 


















X 


V 




0 


















CC 






a 


















3 


c 




0 


















% 


0 






















rH 


*H 




rH 






w 








A 




H 


Z 




rH 






X 








w 






u 




0 






a 








w 




cn 


0 










\ 








'w' 




Ui 


0 




V 






A 








< 




cc 


p 




0 






Ut 








11 




0u 












1 








^-s 




ft 


V 




0 




CM 


rH 








M 




rH 


0 




0 




\ 


w 








M 




> 






rH 


Lu 


< 


* 


A 






v-r 




ft 


z 




0 


* 


X 


X 


CO 






X 




< 


0 




£ 


<c 


* 


o 


>✓ 






• • 




x 


0 






o 


Uu 


cn 


< 






<rft 




cu 


z 


U 


0 


* 


* 


* 


* 




in 


M 




z 




3 


-P 


ui 


Ui 


ui 


kD 




o 


M 




< 


z 


* 


0 


z 


J 


j 


IV 


O H 






ft 


p 


cn 


rH 


o 


O 


o 


• 


• 


rH 


< 






0 


> 


P 


z 


z 


z 


CO 


o 


11 


II 


M 


z 


> 


o 


a 


11 


II 


11 


11 


11 


W 


/A 


M 


M 


c 


it 


rH 


A 






Aft 


A 


W 


M 


H 


CC 


0 


cn 


0 


rH 


(M 


CO 


TJ* 


in cc 


W 


X 


cu 


u 


> 


o 


V 






V 




O 


w 


Ui 


z 


ft 


a 


\ 


< 


< 


< 


< 


< 


Ui 


CQ 


z 


o 


in 


o 


in 


o 


O 


o 


o 


o 


o 


o 


o 


uo 




n 


n 


00 


CT» 


o 


rH 


CM 


CO 




in 


CM 


CM 


CM 


CM 


CM 


CM 


CO 


CO 


CO 


CO 


CO 


CO 



79 



355 'V, II, & 12 are values that are fed into subroutines 5500 & 6000 

356 'These subroutines are used to calculate internal energy, enthalpy and moles 



357 'The temperatures are raised to powers in the polynomial expressions. 

358 'The values of II S. 12 tell the subroutines for what species to solve 
360 Y=TS: 11=1 : 12=4 

370 FOR 1=11 TO 12 
380 GOSUB 5500 
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560 KK= 1 

565 'If injection has occurred, then go to Combustion Subroutine. 
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4360 U<1, 1)=3.0959 
4370 U(l,2)=2. 73114E-03 
4380 U<l,3)=-7. 80542E-O7 
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4750 AC=4»VC/D 
4760 AS=S“PI *D 
4770 AT= AS+AC 
4780 RETURN 



VOLUME & AREA OF CYLINDER 
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5100 ' 

5X01 '********** fuel calculations 

5102 ' ********************************* 
5105 •'Calculate Number moles of fuel. 
5110 S0X=CA+HA/4 
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5200 ' ************************************************** 

5201 '******»»*» PRINT SUBROUTINE ****«***»* 

5202 '*m***»»m****muum**u*****m»m***»****tt*u*****»***»»u 

5205 LPRINT CHRS < 27) ; "6” : LPRINT CHRS <27) ; CHRS (49) ; CHRS ( 30) 

5210 LPRINT:LPRINTTAB<25>"DIESEL ENGINE COMBUSTION CYCLE” : LPRINT 
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5500 '******************************************* 
5505 '****** THERMODYNAMIC PROPERTIES OF MIXTURES 
55X0 'it****************************************** 
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6090 F3=F3+X ( I ) 
6100 NEXT I 
6120 RETURN 
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6670 RRCT=RACT+ <R1N*ADELT> 

6680 IF RRCT < PPEP THEN DF=R1N*ADELT ELSE DF=PN1 * ADELT : GOT06691 



8690 DF=PN1 * ADELT 

6691 DMF=DF/ <WF*M0LE) 

6692 IF (KO+DMF) >1.0 THEN DMF=(1.0-K0> 

6693 IF (KO+DMF) =1.0 THEN ZZ=1 

6700 'Calculate the amount of fuel burned in this step 
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HEAT TRANSFER- ANN AND EQUATION 
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